AD  A 0 5 1 6 5 3 


. • ■ . . 'i-;  VW-.T.  ’ 


• V*;> Yi.“’'>*:^  fA 


Industrial  S Systems 
Engineering  Department 
University  of  Florida 
Gainesville,  FLi  326i  I 


DISTRIBUTION  ST  ■■  T 

Approved  for  prT  T-- 
£>i:tTibuur;'i 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE,  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 


REPRODUCE  LEGIBLY. 


A DYADIC  AGE-REPLACEMENT  POLICY  FOR  A PERIODICALLY 
INSPECTED  EQUIPMENT  ITEM  SUBJECT 
TO  RANDOM  DETERIORATION 

Research  Report  No.  77-6 

by 

Srlnivas  N.  Iyer 
Boghos  D.  Slvazllan 

July  1977 


Department  of  Industrial  and  Systems  Engineering 
University  of  Florida 
Gainesville,  Florida  32611 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED 


I 

’I'l 

f!? 


I-- 


Thls  research  was  supported  In  part  by  the  Army  Research  Office, 
Triangle  Park,  NC,  under  contract  number  DAHC04-75-G-0150. 


THE  FINDINGS  OF  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN  OFFICIAL 
DEPARTMENT  OF  THE  ARMY  POSITION,  UNLESS  SO  DESIGNATED  liY  OTHER 
AUTHORIZED  DOCUMENTS.  ^ 


D D ^ 

pnnjiEi 

22  tsra 


MFE 

B 


i 


H 

-:r 


.. «— 


REPORT  DOCUMENTATION  PAGE 


SBCUBITV  CLAMiriCATlON  OF  THI»  PAOf  flWi**  Bntnmd)  ^ (s  4“^  » ^ ^ 

I" ' "' '-'""'"'r 'i^---*  I ‘ "'r.a#miiu 


HjnnBsTRUCTloNT' 

BEFORE  COUCPLETINO  FORM^ 

mwn,M„wuTim  f.TA. 


A Dyadic  Age-Replacement  Policy  for  a Periodicall;j 
Inspected  Equipment  Item  Subject  to  Random 
Deterioration,  ' ■ — ' 


T.  »UTHQHr»l_ 


Srinivas  N. /lyer 
Boghos  D.ystvazlian 


ir  FCNFONMINO  OROANIIATION  NAME  AND  AbORIII  I..  , 

Industrial  and  Systems  Engineering'/ 
University  of  Florida 
Gainesville,  Florida  32611 


Qg  WO«T  A.Mmi 


Technical 


«.  FKNPONMINO  ONO.  nePONT  NUMBCN 

77-6 


V "--V  ■■  doMTWACT  OH  OWAWT  NUM B « NftJ 

( ( / DAH-C04-  75-0-^15^1 


tt»VBoW«gtV 

zqrf611j(2A14D  /: 


2qB611j(2A14D 
“Ipp'I  of  Applied  Math. 


NUMBBNI 

Rsch  in 


tl.  CONTnOLkINO  OPPiei  NAMI  AND  ADORtIt 

U.  S,  Army  Research  Office 
P.  0.  Box  12211 
Triangle  Park,  NC  27709 
I.  monitobino  aokncy  namb  b add 


U.  MONIToniNO  AOKNCY  NAMB  B AODAMVIf  itihtml  /nm'OmiftHng  Oltlao) 


tit,  pjCLAIIIPICATION/DOWHOPADINa 

KhKulb 


Ur 'OOTRIBuVidN  tVArKMINT  r«r<Al«  Rtpw<> 

Approved  for  public  release;  distribution  unlimited. 


<7.  DIITRtBUTION  ITATtMINT  (•/  Hit  tItIttM  19,  II  mUnml  hm  Rappf^ 

N/A 


iB.  •upplkmkntarv  notes 


t|.  KBV  WORDS  fOanSmM  an  iwvara*  tUt  II  tiMCPMip  anW  Htnlllf  Sr  SI*aS  maiSw; 

Replacement  System 
Aging  and  deterioration 
Linear  Dlophantine  equations 
Economic  replacement  model 
Optimum  dyadic  replacement 


as.  A 


ea^RAeT  (Cmllnt*  m rsrwM  tIH  II  tttttimr  tnd  Umtlllr  Sr  M»oS  mmitt) 


iA  periodic  review  replacement  system  Is  considered.  The  amount  of 
deterioration  over  successive  periods  form  a sequence  of  1.1. d.  random 
variables.  A replacement  policy  of  the  dyadic  type  is  in  effect  whereby  the 
used  equipment  Item  Is  discarded  and  Immediately  replaced  by  a new  identical 
equipment  item  If  at  the  end  of  a period  the  old  equipment  has  service  aged 
by  an  amount  In  excess  of  S or  has  been  In  operation  for  'exactly  N periods 
whichever  comes  first.  Expressions  for  the  joint  distribution  of  the  service 
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operation,  optimal  values  of  S and  N are  computed  for  a few  numerical  examples. 
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Abstract 


A periodic  review  replacement  system  is  considered.  The  amount  of 
deterioration  over  successive  periods  form  a sequence  of  1.1. d.  random  variables. 
A replacement  policy  of  the  dyadic  type  Is  In  effect  whereby  the  used  equipment 
Item  is  discarded  u.\d  lr.u?edlately  replaced  by  a new  Identical  equipment  Item 
if  at  the  end  of  a period  the  old  equipment  has  service  aged  by  an  amount  In 
excess  of  S or  has  been  In  operation  for  exactly  N periods  whichever  comes 
first.  Expressions  for  the  Joint  distribution  of  the  service  age  and  the 
chronological  age  and  for  the  distribution  of  the  total  number  of  replacements 
are  derived.  The  derivation  of  the  distribution  function  of  N|.  relies  on  the 
solution  to  a system  of  linear  Dlophantlnc  equations.  Finally,  using  as 
criterion  the  minimization  of  the  total  steady-state  expected  coht  per  period, 
consisting  of  a fixed  replacement  cost  and  a linear  cost  of  operation,  optimal 
values  of  S and  N are  computed  for  a few  numerical  examples. 
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1.  Introduction 


We  consider  a periodic  review  replacement  system  of  a piece  of  equipment 
item  which  ages  while  in  operation.  The  amount  of  deterioration  or  service 
aging  {D^}  over  successive  equally  spaced  periods  (1  ■>  1,  2,  ...)  form  a 
sequence  of  Independent  and  identically  distributed  positive  random  variables 
with  know  distribution  function  and  p.d.f.  4>(OiO  < ^ We  assume 

that  a replacement  policy  of  the  dyadic  type  la  In  effect  whereby  at  the 
end  of  each  period  the  total  amount  of  deterioration  of  the  equipment  item 
since  acquisition  (service  aging)  la  measured  as  well  as  the  >.0" number  of 
periods  elapsed  since  such  acquisition  (chronological  aging) ; the  used 
equipment  item  is  discarded  and  Immediately  replaced  by  a new  identical  equip- 
ment item  If  at  the  end  of  a period  the  old  equipment  has  service  aged  by  an 
amount  In  excess  of  S or  has  been  In  operation  for  exactly  N periods  whichever 
comes  first.  We  shall  denote  this  policy  aa  an  (S*  N)  pollc/. 

Age  replacement  problems  have  been  studied  by  several  authors  in  the  past 
([!]>  12J|  [4],  [9J|  [11]),  although  these  chiefly  center  around  deterioration 
and/or  breakdown  in  continuous  time.  In  discrete  time  maintenance  models,  the 
approach  has  been  mostly  based  on  Markov  decision  theory  [3],  [6],  [7].  In  this 
respect,  the  model  considered  in  this  paper  departs  from  practice  in  that  It 
studies  a two  decision  variable  (dyadic)  periodic  review  age  replacement  problem 
using  renewal  theory. 

Consider  for  example  an  equipment  item  which  is  operating  intermittently 
characterized  by  the  fact  that  over  equal  time  intervals  the  total  service  provided 
varies  depending  upon  the  user's  request  for  service.  Two  Victors  determine  the 
deterioration  of  the  equipment  Item,  namely,  its  actual  usage  and  Che  time 
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elapsed  since  acquisition.  The  equipment  may  be  a copier  in  a duplicating 
office.  Its  usage  over  successive  intervals  of  months  is  measured  by  the 
total  number  of  reproductions  accumulated  during  the  month.  A maintenance 
policy  is  then  dictated  simultaneously  by  the  total  number  S of  reproductions 
since  last  maintenance  (approximated  by  a continuous  variable)  as  well  as  the 
time  N elapsed  since  such  maintenance.  Another  example  is  provided  by  a road 
vehicle  in  which  the  usage  is  measured  by  the  total  number  of  miles  registered 
during  successive  unit  time  intervals.  The  age  of  the  road  vehicle  is  both  a 
function  of  the  total  number  of  miles  registered  as  well  as  the  time  elapsed 
since  acquisition.  The  road  vehicle  is  replaced  whenever  its  usage  exceeds  a 
given  level  S or  whenever  it  has  been  in  operation  for  a certain  length  of  time 
N whichever  comes  first. 

We  assume  that  at  time  origin,  t 0,  the  equipment  item  has  Just  been 
replaced  and  has  service  and  chronological  ages  equal  to  £.ero.  Following  a 
decision  at  beginning  of  period  t,  t ••  1,  2,  let 

■ total  number  of  replacements  In  the  time  interval  (0,  t] 

■■  service  age  of  equipment  Item,  0 £ Y^  < S 
0^  ■ chronological  age  of  equipment  item, 

Oj.  “ 0.  1 N - 1 

In  what  follows  we  shall  consider  the  two  stochastic  processes 
{(Yt,  Gj.) , t ■ 1,  2,  ...}  and  {N^,  t 1,  2,  ...}  and  derive  in  particular  expressions 
for  the  Joint  distribution  function  of  Y^  and  and  for  the  distribution  function 
of  N^.  Since  the  case  when  N ~ 1 is  trivial,  we  restrict  ourselves  in  the  sequel 
to  the  case  when  N >_  2.  In  deriving  the  distribution  function  of  N^.,  the  analysis 
will  rely  on  the  solution  to  a system  of  linear  Dltphantlne  equations.  Finall:^,  using 
as  criterion  the  minimization  of  the  total  steady  state  expected  cost  per  period,  con- 
sisting of  a fixed  replacement  cost  and  a linear  cost  of  operation,  optimal  values  of 
S and  N are  computed  for  cases  when  {D^}  have  gamma  distributions. 
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2 , Distribution  of  the  Waiting  Time  Till  the  Replacement. 

Let  {Tj},  j “ 1,  2,  ...  represent  the  sequence  of  interarrival  times  between 
replacements.  {Tj } is  an  ordinary  renewal  process  over  discrete  times  [8]. 

The  distribution  function  of  Tj  is  given  by 
P{Tj  ^ 0}  • 0 


P{T^  < n}  - P{D,  + D.  + •..  + D > S} 
j X A n 


1 - , n-l,2,...,N-l 


(1) 


P{Tj  ^ N}  ■ 1 


where  we  define 


- /g  - u)  , n - 1,  a,  - 1) 

Let  f_  (n)  be  the  probability  mass  function  of  the  Interarrival  time  T, . 

J ^ 

It  immediately  follows  from  (1)  that 


(0), 


(2) 


n - 1,  2,  , N - 1 


(3) 


n ■ N 
otherwise 


and 


EIT.]  - I nf  (n) 

J -m 


r-N-1  . , 

1 + I 

r“l 


(4) 


Also,  the  probability  generating  function  of  T, , G_  (u)  , may  be  evaluated! 

^ J 


G„  (u)  - ^ u”  f-  (n) 


T, 

” j 


u < 1 


(5) 


n-N-1  . V 

u-  (1  - u)  I u“  *^"^(8) 
n-1 


^ - 


I - 1 ■ . ■ .» ^ 


Let  Wj^  be  the  waiting  time  till  the  replacement  k ■«  1,  2 The 

probability  mass  function  of  is  clearly 

fy  (t)  - f*  (6) 

k 

where  * denotes  the  usual  convolution  operation.  Thus,  the  probability  gener- 
ating function  of  Wj^,  (u)  is 

(u)  - [u  - (1  - u)  I u*'  k - 1,  2,  ...  (7) 

k n-1 

from  which  the  distribution  function  of  W,  can  be  derived. 

k 

3.  Joint  Distribution  Function  of  Y.  and 
t t 

Let  for  t - 1,  2,  ... 

■>  probability  of  a replacement  being  made  at  beginning  of  time  t. 

0)<iy  " probability  that  at  beginning  of  time  t,  the  service  age  of  the 
equipment  item  lies  between  y and  y -f  dy  and  its  chronological 
ago  is  exactly  0,  0 < y < S and  0 ■ 1,  2,  ...,  N-  1. 

Thus, 

M - P(Y  ■ 0}  0 - 0}  (8) 

t C t 

00  <30  . 

- I “ t}  - I f*  ^*^^(t),  t - 1,  2,  .. . 

k-1  k-1 

Also,  for  0 < y < S and  0-1,  2,  ...,N-1 

'f^(y,  0)dy  - P{y  < 1 y + dy;  0^.  - 9}  (9) 

00 

“ I P{y  < Yj.  ly  + dy|Wj^  - t - 0}.  P {Wj^  - t - 6} 
k-1 

OO 

- I d <J^®^(y)  f!  ^‘'^t-0)  , t - 1,  2,  ... 

k-1 

Some  well  known  results  in  renewal  theory  [8  ] can  be  obtained  from  relations 
(3)  to  (9)  by  considering  the  special  case  when  N equals  to  *. 
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4.  Limiting  Distribution  of  9) 


It  is  easily  verified  that 


M “ lim  M 

t-Ha 


1 

ElT^] 


1 

1 + I 

r-1 


(10) 


To  determine  4'(y,  9)dy  - lim  4*  (y,  0)dy,  where  0 < y < S and 

t-H»  ^ 

6-1,  2,  N - 1,  define  the  generating  function 

OQ 

G(u)  - I (y,  e)dy  |u|  £l 

t-1 

Using  (9)  we  obtain 

G(u)  - I I d f!  ^‘'^t  - 0) 

t-1  k-1 

-u“(d  C-i^gVuT-lI 

Then, 

ny,  0)dy  • lim  (1  - u)u®,[d  ] [-CQ-y-y  - 1] 

u+1 

E[Tj] 

- — , 0 < y < S,  9 - 1,  2 N - 1 (11) 

1 + ^ 
r-1 

Limiting  distribution  under  an  (S,  <«)  policy  follows  immediately  from  (10) 
and  (11).  It  is  also  possible  to  determine  from  (10)  and  (11)  the  marginal  dis- 
tribution of  service  age  at  beginning  of  period  and  end  of  period. 

5.  The  Distribution  of  N^,  t - 1,  2,  ... 


Let  (z,  t)  be  the  probability  generating  function  of  N^,  l.e.. 


5 


(z.  t)  - J P{N  - r}  , Izl  ^ 1 
r-0 


Also,  let  G(z,  w)  be  the  generating  function  of  G (z,  t)  with  respect  to 


t , i.e. , 


G(z,  w)  - ^ G (z,  t)  , |w|  < 1 

t-1  t 


Then,  [8l 


G(z,  w) 


1 - G™  (w) 

J 


(1  - w) [1  - z G (w)  J 


and,  using  (5),  we  obtain 


K-1 


G(z,  w)  - [1  + J ^ I I {z  G_  (w)}  "] 


n"l 


m"0 


5 , 1 Inveraiou  of  G(z.  w) 


In  order  to  determine  P{N^  ■ r}  , r - 0,  1,  2,  . . . , we  first  Invert 


G(z,  w)  with  respect  to  w and  the  resultant  expression  with  respect  to  z. 
For  notational  convenience,  lot 


b^  - n - 1,  2,  ... 


^0  " ^ 


Then  (5)  becomes 


N-1 


L (u)  - [ I - b^)J  + 

J i“l 


and  (15)  becomes 


N-1 


G(z,  w)  " [bg  + I t I z™  { G^  (w))  “] 

n"l  m“0  j 


To  derive  the  coefficient  of  w*^  in  G(z,  w) , t ■ 1,  2,  ...,  it  is  necessary  to 


develop  a systematic  procedure  for  determining  the  coefficient  of  w in  the 


2 N ' 

multinomial  expansion  of  an  expression  of  the  form  + ***  + ) 


where  the  a^,  1 ■ 1,  ...,  N,  are  known  constants,  u is  any  positive  integer, 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 
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and  T = 1,  2,  ...,  uN. 

From  the  Multinomial  Theorem,  the  coefficient  of  is  given  by 

, "l  "2  "n  f 

i-T-f— r-l  “2  "■  “»  ‘ 

^ r>2^!n2!  . . . 

where  the  summation  is  taken  over  all  non-negative  integers,  n^,  ± ■ 1,  2,  N 


such  that 


+ n2  + • • • + Hjj  “ u 


n,  + 2n„  + •••  + Nn„  ■ t (2 

12  N 

« 

If  no  set  of  n^,  i ■ 1,  2,  N exists  satisfying  (20),  (21)  then  the 

coefficient  is  equal  to  zero. 

5.2  On  the  Solution  of  a System  of  Two  Linear  Diophantine  Equations 

(20) , (21)  represent  a system  of  two  linear  Diophantine  equations  in  N 
non-negative  Integral  unknowns.  It  has  been  shown,  [5],  that  any  system  of 
m linear  Diophantine  equations  can  be  transformed  into  a single  aquation  and  a 
set  of  inequalities.  For  ovr  special  structure,  the  following  interesting  theorem 
is  obtained  by  applying  the  so-called  'Rule  of  Virgins'  (method  of  eliminating 
one  unknown  at  a time)  [ICj . 

Theorem;  The  system  of  two  linear  Diophantine  equations  in  N ( > 1)  non-negative 
integral  unknowns 

n,  + n-  + • • • + n - u (2 

12  N 

n^  + 2n2  + *••  + Nn^j  - t (5 

where  u and  x are  positive  integers,  has  a solution  if  and  only  if  u <_  i u N. 

If  it  has  a solution,  then  all  possible  solutions  of  (20),  (21)  are  generated 
by  all  possible  solutions  of 


I 


c + c-o  + •••  + c u - T + (N  - 1)  [t/N]  (22) 

1 Z N“  i 

and 

[t/N]  - u ^ £ Cj^_2  1 • • * i.  *=2  ~ 

where  the  i“li  2,  ....  N - 1 are  integers  (unrestricted  in  sign)  and 
[t/N]  is  the  greatest  integer  £ t/N;  and  the  n^,  i • 1,  2,  ....  N,  are 
given  by 

n^  - 0 - [t/N]  + (24) 

njj_  j " ” °J+1  ...,N-2 

- [t/N]  - 

Further,  If  the  original  system  of  equations  has  no  solution,  then  there  is 
no  solution  to  (22),  (23),  so  that  the  two  systems  are  equivalent. 

Proof ; It  is  easily  verified  that  if  t < u or  if  t > u N,  then  there  is  no 
solution  to  the  system  of  two  linear  Dlophantlne  equations.  Also,  if 
u ^ T _<  y N,  it  is  clear  that  the  coefficient  of  w^  In  (a^^w  + + •••  + *j^**)*^ 

is  given  by  (19)  (this  coefficient  being  not  necessarily  non-*ero) ; so  that 
there  is  at  least  one  solution  to  (20),  (21). 

To  prove  the  second  part  of  the  theorem,  assume  that  there  exist 
n^,  i-  1,  2,  ...,  N satisfying  (20),  (21).  Then  from  (21)  we  have, 

N-1 

n„  - t/N  - ( ^ in.)/N 
” i-1  ^ 

N-1 

■ [t/N]  -(.  I in.  + N[t/N]  - t)/N 
i-1 


Now  let 

N-1 

c - ( ^ in.  + N[t/N]  - t)/N 
^ i-1  ^ 

Clearly,  c^  is  an  integer,  and  n^^  being  non-negative,  £ [t/N].  Also, 
Hj,  - [T/N]  - C,^ 


(25) 


(26) 


1 


! 
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Again,  from  (25),  we  have. 


N-2 

”n-1  " ~ ^ ^ “ ■r)/(N-l) 

N-2 

Letting  S " ^ I S N[t/N]  - t)/(N-1),  we  get 

^ i-1 

Vi  “ *^1  " 

Clearly,  is  an  integer,  j< 

Proceeding  in  this  fashion,  we  obtain 

"(N-J)  " " ‘^j+l  J 2 n-2 

N-j  J-1 

where  c “ ( |[  In^  " I S N(t/N)  - t)/(N-J  +1)  J-2,  3,  ...,N-1 

J 1-1  ^ 1-1 

and  Cj  ^ 1 ^ ••••  N — 1 

Letting  ,j  - N - 1 In  (28) , we  obtain 
N-1 

n^  - - N[t/N1  + T 

From  (20),  (26),  (27)  and  (29),  we  have, 

N-1 

I c - u - T + (N  ~ 1)  [t/N] 
i-1  ^ 

so  that,  from  (29), 

nj^  - u - [t/N]  + 

and  since  n^^  ^ 0,  [t/N]  - 0 £ “^N-l 

Thus,  starting  from  (20),  (21)  we  have  obtained  (22),  (23)  and  (24).  Finally, 
it  can  be  easily  shown  that  any  solution  to  (22),  (23)  yields  a solution  to 
(20),  (21)  given  by  (24),  so  that  the  two  systems  are  equivalent. 


To  illustrate  the  use  of  the  theorem  we  solve 


01  + 02  + 03-5 

+ 202  ^"3  " ® 

Here,  N-3,  u-5,  t-8  and  [t/Nl  - (8/3 J - 2 
From  (22),  (23), 

Cl  + C2  ■ 5 - 2 - 2 - 1 
subject  to  £ Cg  £ < 2 

which  yields  as  solutions,  (c^  - 2,  C2  - -1}  aod  {cj^  " li  C2  - 0}. 
The  corresponding  solutions  to  the  original  system  are  from  (24) : 


{n^^  - 2,  02  ■ 3,  Oj  - 0}  and  {n^^  - 3,  02  ■ 1,  Oj  » 1} 


It  should  be  noted  here  that  the  computational  effort  required  to  solve 
the  transformed  system  is  no  greater  than  that  for  the  original  system.  In 
fact,  In  many  Instances,  especially  when  the  nuid>er  of  solutions  is  not  large, 
the  transformed  system  is  easier  to  solve,  being  more  amenable  to  an  snunsrative 
scheme.  Thus,  although  the  theorem  presented  is  interesting  In  itself,  It  also 
possesses  some  utility. 

5.3  An  Expression  for  P{N^  ■ r}  t ■ 1,  2,  ...  ; r - 0,  1,  2,  ... 


Let  I be  the  number  of  solutions  to  the  system  of  linear  Diophantine 


U,T 


.th 


equations  (20),  (21).  Let  {n^^  , 1 ■ 1,  2,  . . . , N}  be  the  J solution, 


j - 1,  2,  ...,  if 


> 0.  Define 


„ "2j  . ^ "nj 

a.  ‘'a„  •'  •••  a^ 

1 4 n 


u.T,j 


"lj‘  "2j' 


"nj 


r * J " 2 


> 0 


(30) 
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X <2  N O 

Then,  from  (19),  the  coefficient  of  w In  (a^^w  + BjW  + •••  ) Is  given  by 


(■I'l 

m 

‘ ;■! 

d 

'.‘■I 


? f ■ 


['  ’ 

c. 

J 

\ 

h 

k' 


i 

» 

&■ 


lu  ,T 


I A , 
j-1 


If  I ■■0,  the  above  expression  ia  defined  to  be  sero. 


U.T 


Ut 

- b^_j  " **1  • i ■ 1.  2,  . .. , M - I 

*N  ■ Vl 

Then  (17)  becomes 

and  hence  it  follows  from  (31)  that  the  coefficient  of  , t ••  1,  2,  . . . in 

[G«  (w)]*',  k - 1,  2 is  given  by 

J 

'k.t 

^ t J 

j-1  ’ 

From  (18),  it  is  easily  seen  that  when  1 < N £ t,  the  coefficient  of  w^  in 
G(z,  w)  is  given  by 

'I  b f {Coefficient  of  w**  in  [G_  (w)  ]*^} 
q-t-N+1  k-1  ^j 

and  for  all  values  of  N greater  than  t,  the  coefficient  of  is  given  by 

I {Coefficient  of  in  [G_  (w) 

^ q-1  ‘•“‘I  k-1  j 

Thus,  from  (34),  (35),  (36),  we  get 


I PCN 
r-0 


t g . k,q 

- r)  - I b Ikl*  I A.  .l<N<t 
q-t-N+1  ^ **  k-1  J-1 


and 
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(31) 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


N > t 


(38) 


I P{N  - r}  - b + I b ? kl  t \ q j 

r-0  ‘ q-1  *=  k-1  J-1  ^’‘*•3* 

To  derive  an  expresalon  for  P{N^  '>  r),  t *•  1,  2,  ...,  r ■ 0,  1,  2,  ...,  wc 

will  distinguish  between  two  cases,  one  where  1 < N £ t and  the  other  where 

N > t 

Case  I 1 < N < t 


From  (37),  we  have 


P{N^  - 0}  - 0 


P{N^  ■ r} 


P{N^  - r}  - r! 


(39) 


, r ■ t 1,  t + 2,  ... 

S V*! 

q-t-N+r  •="‘1  J-1 


IT  - 1,  2, 


Case  2 N > t 


From  (38) , ve  have 


P{N^  - 0}  - b^ 

P{N^  -r}-0  ,r-t+l,  t+2,  ... 


(40) 


t r,q 

P{N.  - r}  - r!  I b_  I A 
t qir  j-1 


1.  2, 


It  Is  clear  from  (40)  that  the  distribution  of  the  number  of  replacements 
In  a time  Interval  (0,  t]  Is  the  same  for  all  values  of  N greater  than  t.  In 
this  case,  while  solving  the  Dlophantlne  equations  (20),  (21)  with  r and  q as 
parameters  Instead  of  u and  t,  we  can  conveniently  set  N ■ t -f  1 in  order  to  have 
a mlnlmvim  number  of  unknowns,  although  any  value  of  N > t would  yield  tha  sans 
final  probability  given  by  (40) . 

5.4  A Sample  Calculation  For  t - 10,  N - 3 

Note!  b as  usual  denotes  «^"^(S),  and  tha  results  are  trua  for  any  distribution 
n 

^(O»0<C<"  and  any  S 
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■ r}  - 0 , r - 0,  1,  2 


r}  - 0 , r ■ 1.1 , 12,  ... 


3}  - - bj)  + bj^b^ 

A}  - bj,(6<l-b^)^  b^  + 12(1  - \)0>i  - bj)^  bj  + (b^  - bj)^! 

+ b^(12(l  - b^Xb^  - b^)  b^  + U(h^  - bj)^  bj]  + 

[4(1  - b^)  hi  + 6(bj  - bj] 

5}  - b2[20(l  - b^)^  (b^  - bj)  bj  ■*•  10(1  - bj>^  (b^  - bj)^]  f 
+ b,  (10(1  - b^^  bl  + 30(1  - b,)^  (b,  - b,)^  b,  + 


1/  - ^2'  “2 


5(1  - bj)(b^  - bj)^]  + [30(1  - (b^  - bj)  b|  + 

20(1  - b^Xbj^  - bj)^  bj  + (bj  - bj)^] 

6}  - b2[6(l  - bp^  bj  + 15(1  - bj^)^  (bj^  - bj)^) 

b^[30(l  - bj)^  (b^  - b2)b2  + 20(1  - b^>^  (b^  - bj)^!  + 


'1'  '•'1  “2' 


[15(1  - bp^  b2  + 60(1  - b^>^  (bj  - 62)^  62  + 

15(1  - b^)^  (b^  - b2)^] 

7)  - b2(7(l  - b^)^  (b^  - b2)]  + b^[7(l  - b^)*  b2 

21(1  - h^)^  (b^  - b2)^]  + [42(1-  b^)’  (b^  - bj)  b2  + 

35(1  - b^)^  (bj  - b2)^] 

8}  - b2(l  - b^l®  + b^(8(l  - b^)^  (b^  - b2)]  + 

[8(1  - bj^)^  b2  + 48(1  - bj)®  (b^  - b2)^l 
9}  - bj^[(l  - bj^)®]  + [9(1  - bj)®  (bj^  - b2)l 


P{Nj.  - 10}  - (1  - 


I; 


5.5  Tlie  Mean  Value  of  the  Number  of  Replacemanta  In  a time  Interval  (0,  t] 
The  formal  expressions  obtained  from  (39)  are 

Casa  1:  M < t 


B(N^]  ■ I rP{N  - r}  - J r^Cr  - 1)J  J 

“ . ^ . 


r«l 


r-1 


q“t-N+l 


t-q 


I A, 


J-1 


r.q.J 


(^1 


Casa  2:  N > t 


E[N 


] - I rP{N  - r}  - I r-'Cr  - 1)!  I I A 

r-1  ' r-1  q-r  j-1 


j 


(42 


6.  An  Economic  Raplscanant  Modal 

In  ehls  sactlon  wa  consider  tha  problan  of  minimising  tha  total 
staady-atata  axpactad  cost  per  parlod,  which  Is  the  sum  of  the  steady-state 
axpactad  costs  par  period  of  raplscanant  and  operation.  In  formulating 
this  cost  objactlva  functlon»  S and  N appear  as  tha  decision  variables  and 
optinal  values  of  thaao  will  ba  obtained.  It  will  be  assumed  that  wa  have 
a fixed  recurring  cost  of  replacauant  per  period  and  a linear  expected  oper- 
ating cost  par  period.  The  analyals  will  be  made  for  the  case  where  the 
service  aging  par  period,  (i  - 1,  2,  ...},  has  a gamma  distribution.  That 
for  the  negative  exponential  distribution  will  be  derived  as  a special  case. 

To  this  end,  we  Introduce  tha  following  notation: 

LCy,  9)  * conditional  expected  operating  cost  per  period  given  chat  at 

Che  beginning  of  each  period  the  equipment  item  has  chronologlal 
age  e,  (6  - 1,  2,  N - 1)  and  service  age  y (0  < y ^ S) 

L(0,  0)  - conditional  expected  operating  cost  per  period  given  that  the 
equipment  item  has  just  been  replaced. 

K - fixed  cost  of  replacement 
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R(S,  M)  * ■taady^atats  expected  cost  of  replacement  per  period 
0(S,  N)  • steady-atate  expected  coat  of  operation  per  period 
F(S,  N)  ■ ateady-atate  total  expected  coat  per  period 
- R(S,  N)  + 0(S,  N) 


The  reaulta  (10)  and  (11)  obtained  In  Section  4 for  the  ateady-atate 
probability  of  a replacement  being  made,  M,  and  for  the  limiting  dlatrlbutlon 
of  {(Y^,  0^):  t ■ 1|  2,  ...}  will  now  be  uaed  hare  in  formulating  the  coat 
objective  function,  F(S,  N). 

We  have,  from  (10) 


R(S,  N)  - KM  •f  ML(0,  0) 


« K ^ L(0.  0) 
0-N-l 

I *(«) 

e-]  • 


(S) 


(43) 


and  from  (11) 

0(8,  N) 


0-M-l  f® 

I _ L(y,  0)  Y(y,  0)dy 

0-1  Jo 


0-U-l  r 

I L(y,  0)  /^®\y)dy 

0-1  Jo 


0-N-l  .-V 
I + I 0^®^S) 
0-1 


(44) 


where  ^'*^®^(y)  denotea  ea  uaual  the  6-fold  convolution  of  the  denalty 
function  i^(y)  correapondlng  to  the  dlatrlbutlon  function  t(y). 

From  (43)  and  (44),  we  have. 


F(S,  N) 


K + L(0,  0)  X 
8-1 


9-N-l  ( 


Jo 


L(y,  0)  /^®^y)dy 


‘ -0  * O”-” 


(45) 
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AatunMi  a llnaar  conditional  expected  operating  cost  in  y and  e,  and  let 


L<yi  6)  - A + By  + C0  , 0 < y £ S;  6 - 1 , 2 , . . . , N - 1 
L(0,  0)  - A 


Let  the  eervlce  aging  per  period  have  a ganraa  distribution  with 


parameters  r and  1 l.e.  let 


♦(y) 


y i 0 


0 < y < • 


Then. 


rTO 
^0 


(y) 


y 1 0 


_CXy)_^^^_Xe“^y  0 < y < 

rTre) 


-t  a-l 

Defining  y(»,  x)  to  be  the  incomplete  gamma  function  e t d 

Jo 


ve  have» 


Jo  ^ (y)‘*y  • r^)  Y<t0. 


and,  from  (46),  after  some  manipulations. 


fS 


L(y,  9)  d*^®Ny)dy  - (A  + By  + C0)  d"''"(y)dy 

^0  Jo 


fS 


^ I(A  J + Ce)  Y(r9,  XS)  - f 


From  (49),  (47),  (50)  and  (51),  we  heve. 


0-N-l 


i 77^  + X Y(r0,  XS)  - 


B (xs; 

X 


When  the  parameter  r is  an  Integer,  then 


. y (xs)^ 

11 


P(re;  XS) 


(53) 


where  P(r0j  AS)  Is  defined  to  be  the  complementary  cumulative 
Poisson  vlth  parameter  A S. 

In  the  special  case  when  r ■ 1,  we  have  a negative  exponential  distri- 
bution, and  (52)  becomes,  after  some  algebra, 

2 

F(S,  N)  - {K  A + AS<A  + C)  + (-  + C)  - 

M A 

2 

[ASA  + ^"1^-  (“  + 0]  P(N  - 2;  AS)  + [(M  - 1)  A - ASC]  P(N  - 1;  AS) 

+[(i  + c)  P(N{  AS)}  / {1  + AS  - AS  P(N  - 2;  AS) 

+ (N  - 1)  P(N  - 1}  AS)  , N > 2 
6,1  Computational  Results 

Owing  to  the  complexity  of  axpreaslons  (52)  and  (54)  for  the  total  steady- 
state  expected  cost  per  period  it  is  analytically  difficult  to  obtain  optimal 
values  of  N and  S that  mlnlmixe  F(S,  N) , this  difficulty  being  further  com- 
pounded by  the  fact  that  one  decision  variable  Is  discrete  and  the  other 
continuous.  In  what  follows  several  computational  results  are  presented  and 
discussed. 

In  Figures  1 and  2 the  total  steady  state  expected  cost  per  period  Is 
plotted  against  discrete  N for  various  fixed  values  of  S (oi;  alternatively,  AS). 
The  resulting  graphs  have  however  been  made  continuous  for  convenience.  The 
parameter  r of  the  gamma  distribution  (48)  Is  taken  to  be  1 and  3,  respectively, 
so  that  Figure  1 actually  represents  the  case  when  the  service  aging  per  period 
has  s negative  exponential  distribution.  In  both  graphs  the  values  of  the 


other  Input  parameters  are  as  follows:  A - 10  , K ■ 5000,  A - 51,  1. 

and  C ■ 300. 

For  the  negative  exponential  distribution  the  optimal  values  of  N and  AS 
are  found  to  be  6 and  17.0,  respectively,  whereas  for  the  case,  r 3, 

(Fig.  2),  these  values  are  6 and  35.0. 

In  general  the  patterns  in  Figures  1 and  2 exhibit  a decreasing  tendency, 
followed  by  an  Increasing  one,  this  in  turn  giving  way  to  a level  stretch. 

For  fixed  AS,  as  the  limit  on  the  chronological  age  is  raised,  the  equipment 
Item  tends  to  be  replaced  after  a longer  period  of  time  and  hence  the  expected 
coat  of  replacement  decreases  while  the  expected  cost  of  operation  increases 
with  N.  Once  N la  raised  beyond  a certain  "high"  value  (this  "high"  being 
relative  to  the  value  of  AS) , however,  It  no  longer  dictates  the  replacement 
policy.  It  Is  the  "lower"  barrier  of  S that  does  bo.  Therefore,  Increasing 
N does  not  affect  the  expected  costs  of  replacement  and  operation  beyond  this 
point.  The  behaviour  of  the  graphs  is  now  readily  understandable  from  the 
fact  that  the  total  steady-state  expected  cost  is  nothing  but  the  summation 
of  the  expected  costs  of  replacement  and  operation. 

For  very  low  values  of  AS,  N Is  irrelevant  to  the  doclslon-maklng  and  the 
level  stretch  dominates  early  In  the  game.  As  AS  is  raised  further,  the  increas- 
ing and  decreasing  portions  of  the  graph  both  begin  to  make  themselves  felt. 

For  very  high  values  of  AS,  the  curves  converge,  because  now  AS  becomes  Irrele- 
vant and  N dominates  the  decision-muklng. 

When  r is  Increased  from  1 in  Flg^  1 to  3 In  Fig.  2 the  optimal  values 
of  N remains  unchanged  at  6 whereas  that  of  AS  Increases  from  17.0  to  35.0. 

This  Is  explained  by  the  fact  that  Increasing  r increases  the  expected  value 
of  the  service  age,  given  a certain  chronological  age  E)  at  the  beginning  of  a 
period. 

In  both  cases  when  N is  kept  at  Its  optimal  value  of  6,  and  AS  is  increased 
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beyond  Its  own  optimal  value  of  17J0  or  35.0,  the  total  steady-state  expected 
cost  per  period  remains  unchanged  because  AS  becomes  "high"  relative  to  this 
value  of  N and  N dictates  the  replacement  policy.  Hence,  it  would  be  tech- 
nically correct  to  say  of  Fig.  1,  for  example,  that  the  optimum  value  of  AS 
is  greater  than  or  equal  to  17.0. 

Both  Figures  1 and  2 were  plotted  against  N for  various  values  of  AS. 

It  is  easy  to  see,  however,  that  the  behaviour  of  the  graphs  would  be  essen- 
tially the  sane  If  the  plots  were  Instead  made  against  AS  for  fixed  values  of 
N. 

Varying  the  other  Input  parameters  would  naturally  affect  the  optimum 
values  of  N and  S.  If  the  value  of  C is  lowered,  for  Instance,  to  30,  (Fig.  3) 
' the  dacraaalng  portion  of  the  graph  dominates,  and  the  optimal  values  of  N and 
' AS  are  found  to  be  as  high  as  18  and  30*0,  respectively.  It  is  obviously 
more  economical  to  run  the  equipment  for  a longer  period  of  time  rather  than 
Incur  a relatively  exorbitant  cost  of  replacement. 
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